The Astrophysical Journal, submitted 

Preprint typeset using L^TgX style emulateapj v. 1 1/12/01 



THE ORIGIN OF POLAR RING GALAXIES: EVIDENCE FOR GALAXY FORMATION BY COLD 

ACCRETION 

Andrea V. Maccio 1 , Ben Moore' & Joachim Stadel 1 

The Astrophysical Journal, submitted 

ABSTRACT 

Polar ring galaxies are flattened stellar systems with an extended ring of gas and stars rotating in a plane almost 
perpendicular to the central galaxy. We show that their formation can occur naturally in a hierarchical universe 
where most low mass galaxies are assembled through the accretion of cold gas infalling along megaparsec scale 
filamentary structures. Within a large cosmological hydrodynamical simulation we find a system that closely 
resembles the classic polar ring galaxy NGC 4650A. How galaxies acquire their gas is a major uncertainty in 
models of galaxy formation and recent theoretical work has argued that cold accretion plays a major role. This 
idea is supported by our numerical simulations and the fact that polar ring galaxies are typically low mass systems. 

Subject headings: cosmology: theory — galaxies: formation — methods: numerical 



1. INTRODUCTION 

Polar ring galaxies are in general SO galaxies that are sur- 
rounded by a perpendicular "polar" ring or disk of stars and 
gas (Schweizer et al. 1983; Sparke 1986; Whitmore et al. 1987, 
1990; Rubin 1994). Less than one percent of SO's have polar 
rings, however these features are typically only seen edge on 
since they are very low surface brightness, therefore the true 
fraction of galaxies with such features is likely to be higher 
(Figure 2(a) shows a classic example). The currently favored 
scenario for their formation is a merger event between two galax- 
. ies, such as the accretion of a large gas rich satellite which is 
disrupted in a ring surrounding the central galaxy (Schweizer et 
al. 1983; Reshetnikov 1997; Bekki 1998; Tremaine & Yu 2000; 
Sparke & Cox 2000; Bournaud & Combes 2003). However, 
the observed mass of baryons in the ring can be comparable to 
the mass of baryons in the central disk implying that a major 
merger must have occurred. Similar mass mergers lead to co- 
alescence of the two haloes and central galaxies, the remnant 
violently relaxating into a spheroidal system elongated in the 
plane of the merger. This is now the standard model for the for- 
mation of elliptical galaxies (Toomre & Toomre 1972; Gerhard 
1981; Barnes & Hernquist 1996; Burkert & Naab 2003) there- 
fore the initial conditions must be fine tuned in order to create 
a polar ring galaxy from a merger event. 

Galaxy formation is a major unsolved problem that brings to- 
gether many aspects of astrophysics. Observations across many 
scales and redshifts, from the cosmic microwave background 
to the large scale clustering of galaxies has defined a standard 
cosmological model which sets the background cosmology and 
initial conditions within which we can try to understand the for- 
mation of cosmic structures (e.g. Tegmark et al. 2004). Within 
this framework, it is well understood how the dark matter struc- 
tures collapse and form galactic haloes, but much more com- 
plex is the behaviour of the baryons. Even well posed questions 
remain unanswered, such as how the gas settles into the centres 
of dark matter haloes (Hoyle 1953; Birnboim & Dekel 2003; 
Katz et al. 2003; Keres et al. 2005). 

The models that have been proposed for how galaxies get 
their gas are quite diverse. The oldest idea is that the baryons 
are initially shock heated during an early merging and virializa- 
tion phase of the dark matter haloes. The gas begins to settle 
quietly into the central halo in a cooling flow as the ionised 
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gas in galactic haloes radiates photons, losing energy and gain- 
ing angular momentum. Observational support for hot gaseous 
haloes around galaxies is lacking (Benson et al. 2003) but see 
(Kaufmann et al. 2005; Mo et al. 2005). More recently, a wealth 
of observational evidence suggests that a significant reservoir 
of cold gas exists within galactic haloes and groups of galax- 
ies (e.g. Westmeier et al. 2005). This cold material may be 
stripped from infalling satellites or be infalling directly via a 
smooth cold flow focused along filamentary structures. 

2. NUMERICAL SIMULATIONS 

Numerical simulations provide the machinery to evolve the 
initial conditions through the non-linear regime. We imprint 
linear fluctuations on the dark matter and baryons using a two 
component particle distribution according to expectations from 
the standard ACDM model. The simulations were performed 
with GASOLINE, a multi-stepping, parallel TreeSPH A/-body 
code (Stadel 2001, Wadsley et al. 2004) which can follow the 
shock heating and radiative cooling processes. We include ra- 
diative and Compton cooling for a primordial mixture of hy- 
drogen and helium. The gas can form stars if the physical 
conditions are appropriate and stellar evolution can feedback 
energy into the interstellar and intergalactic medium through 
supernovae explosions. The star formation algorithm is based 
on a Jeans instability criteria (Katz 1992), where gas particles 
in dense, unstable regions and in convergent flows create star 
particles at a rate proportional to the local dynamical time (see 
also Governato et al 2004). The star formation efficiency was 
set to 0.1, but in the adopted scheme its precise value has only 
a minor effect on the star formation rate (Katz 1992). The code 
also includes supernova feedback as described by (Katz 1992), 
and a UV background following Haardt & Madau (1996). 

We first perform a dark matter only simulation of a cubic re- 
gion of a concordance LCDM universe (A=0.7,£7q=0.3, <78=0.9) 
that is 90 Mpc on a side. From this we selected a smaller vol- 
ume that by a redshift z = forms a large candidate Galactic 
mass halo (M,/„, k, 10 12 M q ) and several lower mass satellites 
and surrounding galaxies. This region is traced back to the 
initial conditions and repopulated with many more lower mass 
particles (c.f. Katz & White 1993), including a gaseous com- 
ponent throughout the entire high resolution region. The mass 
per particle of the dark matter and gaseous particles are respec- 
tively m d = 4.90 x 10 5 M Q and m g = 9.75 x 10 4 M Q . The dark 
matter has a spline gravitational softening length of 100 pc and 
we have about 4 x 10 6 particles for each component (dark and 
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FIG. 1 . — The distribution of gas and dark matter in a 2 Mpc wide slice of our simulation centered on a galactic mass dark matter halo. The gas temperature is 
coloured purple-white in a logarithmic temperature scale from 10 4 — lO 6 ^. The narrow filament of gas flowing into the halo is at a temperature of ~ 2 X 10 4 X\ The 
filament in the dark matter is significantly wider than the gas and consists of a smooth component and many smaller haloes. The entire simulation volume is 90 
Mpc 3 and is centred on a region forming several galaxies which are simulated at a higher resolution than the rest of the cube. The masses and softening of the dark 
matter and gas particles are 4.90 X 10 5 M Q , 9.75 X 10 4 M Q and 0.2, 0.15 kpc respectively. The polar ring would appear face-on in this projection, perpendicular to 
the filament 



gas) in the high resolution region. 

3. RESULTS AND DISCUSSION 

Figure 1 shows the final distribution of gas and dark mat- 
ter within a 2 Mpc cube surrounding the central halo that hosts 
our candidate polar ring galaxy. As pointed out by previous 
authors (Birnboim & Dekel 2003; Katz et al. 2003; Weinberg 
et al. 2002), we find that the dominant mode of galaxy forma- 
tion is via the accretion of cold gas along filamentary structures. 
The gas in the filaments is not completely cold, but is heated to 
15,000K from the collapse and formation of these structures. 
Some of the gas is shock heated to the halo virial temperature 
but most is accreted in this warm phase. In our simulation vol- 
ume with a dozen high resolution galaxies, this is the only sys- 
tem with a long lived polar ring, consistent with the observed 
abundance of these galaxies. 

Figure 2(a) shows the hubble space telescope image of NGC 
4650A, a classic example of a polar ring galaxy (Sersic 1967). 
Recent observations of its kinematics confirms that the ring of 
stars and gas rotates around the central SO galaxy (Arnaboldi et 
al. 1997; Gallagher et al. 2002). Figure 2(b-d) show close up 
images of the galaxy that is forming at the centre of the halo 
in Figure 1 . The central galaxy is a rotationally supported disk 
that could be classified as an SO due to its featureless "thick 
disk" morphology. The galaxy has an extended ring of gas (and 
stars) in near circular motion perpendicular to the disk similar 
to NGC 4650A. The face on view of the ring clearly shows that 
the morphology of the ring is clumpy and diffuse. 

The length of the SO major axis is 5.1 kpc, and the ma- 
jor/minor axis ration is « 9.5. The SO baryonic mass is 2.95 x 
1O 9 M , with 17% of gas. The ring is mainly made by gas, 
the total mass (gas+stars, not considering the SO galaxy) is 



1.23 x 10 9 M Q (inside a cylinder with a radius of 18 kpc and 
a high of 7 kpc, centered on the SO galaxy), with a M g /M s ratio 
of about 5.6. In the same volume the dark matter mass is 6.75 x 
1O 9 M . The dark matter halo of this galaxy is quite extended, 
and its virial radius (defined according to a spherical overden- 
sity A = 95) is 98 kpc, for a total mass of 5.25 x 10 10 M o 

To test test the stability of the ring we cut a box of 9 Mpc 
side around our galaxy and then we let it evolve in the future. 
The ring is still present after « 1 Gyrs, then it is destroyed by 
the dynamical effects of a massive (1.1M Q ) merging satellite, 
if we remove all the accreting satellites around our galaxy the 
PR survives up to 1.6 Gyrs. 

The ring forms from cold gas that flows along the extended 
~1 Mpc filament into the virialised dark matter halo. The gas 
streams into the centre of the halo on an orbit that is slightly 
offset from purely radial infall and is "braked" by ram pres- 
sure drag from the ionised halo gas. As it reaches the centre 
it impacts with gas in the halo of the existing galaxy and with 
the warm gas flowing along the opposite filament. The small 
impact parameter is sufficient to leave the gas with sufficient 
angular momentum that it can rotate at the halo circular veloc- 
ity at 10-20 kpc. The ring is clumpy and where it has a high 
enough density it can fragment and form stars according to our 
prescription of star formation (Katz 1992). 

The kinematics of the SO and the inclined ring can be used to 
measure the mass distribution of the stars, gas and dark matter 
within two planes and thus the shape of the potential can also 
be constrained (Schweizer et al. 1983; Reshetnikov & Combes 
1994; Arnaboldi & Sparke 1994; Sackett et al. 1994; Combes & 
Arnaboldi 1996; Iodice et al. 2003). Figure 5 shows the shape 
of the dark matter halo and stellar component of the simulated 
galaxy. The effect of dissipation alters the shape of the dark 
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matter distribution (Dubinski 1994; Kazantzidis et al. 2004). 
At the radius of the ring the halo is closer to spherical, but the 
potential (from the baryons and dark matter) is oblate since the 
stellar component is rotationally flattened. 

The potential at ^20 kpc is still oblate with a flattening q = 
0.9. Therefore any gas falling away from the perpendicular 
symmetry axis will rapidly precess into the plane of the disk 
(Weinberg et al. 2002; Murali et al. 2002). Only gas accreting 
perpendicular to the major axis of the oblate potential will sur- 
vive for more than a few dynamical times. Indeed, the ring 
is long lived in our simulation. This mode of gas accretion 
will usually be off-axis and provides a source of material from 
which the disk can grow and form new stars. Gas accreting onto 
a bulge or spheroidal would likely produce an SO galaxy. Calcu- 
lations of larger volumes are required to study the statistics and 
properties of polar ring and normal galaxies that form via cold 
accretion. We expect that there is a wide variety in the proper- 
ties of such systems, from the mass ratios to their kinematical 
states. Accretion of cold gas could also explain other strange 
kinematical features in galaxies, such as counter-rotating stel- 
lar and gaseous disks (Merrifield & Kuijken 1994; Bertola et al. 
1996). 

We thank Linda Sparke for motivating us to investigate the 
origin of these fascinating galaxies. We also acknowledge use- 
ful discussions with M. Arnaboldi, J. Bland-Hawthorn, S. Ka- 
zantzidis and C. Mastropietro. James Wadsley has carried out 
most of the development of the GASOLINE SPH code used and 
the simulations were carried out on the zBox supercomputer at 
the University of Zurich. 
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FIG. 2. — Panel (a) The Hubble Space Telescope image of the classic polar 
ring galaxy NGC 4650A (side-length: 24 kpc), courtesy of the Hubble Her- 
itage site, STSCI. Panels (b-d) show the smoothed images of the stars and gas 
of the simulated galaxy in different projections (side-length 35 kpc): (b) in- 
clined 20 degrees (c) the ring is edge on and (d) face on (this Figure has the 
same axis orientation of Figure 1). The face on image shows a clumpy distri- 
bution of gas that lies in two slightly offset rings. 
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FIG. 3. — Rotation of the stars in the polar ring. Face-on projection, same 
of figure 2d but rotated of Ri 75 degrees clockwise, in order to have the galaxy 
parallel to the the y-axis; side-length 40 kpc. The central region is not plotted 
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FIG. 4. — Line of sight velocity of gas (open circles) and stars (solid squares). 
In both cases the slit is along the major axis of the galaxy. In the upper panel 
the ring is seen edge-on, in the lower face-on. In the lower panel is clearly 
visible the spike in the los velocity due to the SO galaxy. 
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FIG. 5. — The shape of the dark matter and gaseous components of the 
simulated polar ring galaxy. The dashed curve shows the shape of the dark 
matter halo in a simulation where the gas is not allowed to cool. The solid and 
dotted curves show the shape of the dark matter and stellar components in the 
full hydro-dynamical simulation. The dark halo becomes rounder but its shape 
stays quite triaxial in the vicinity of the disk. The stars are highly flattened 
indicating the central disk component. 



